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ABSTRACT 

Strain  energy  expressions  are  obtained  for  rings  of 
rectangular,  T-  and  I-  section.  The  expressions  are  intended 
for  use  in  the  dynamic  analysis  of  ring  stiffened  cylindrical 
shells.  The  approach  is  essentially  a generalization  of  the 
conventional , approximate  analysis  of  straight  beams,  i.e. 
the  influence  of  shear  stresses,  and  of  direct  stresses  at 
right  angles  to  the  axis  of  the  beams  is  neglected. 


1. 


INTRODUCTION 


The  results  obtained  in  this  study  will  facilitate  the 
determination  of  the  modes  of  free  vibation  and/or  of  the 
dynamic  response  of  ring-stiffened  cylindrical  shells. 

Expressions  for  the  potential  energy  of  the  stiffeners  in 
their  most  general  state  of  displacement  will  be  derived. 

The  strain  energies  are  obtained  on  the  basis  of  "beam  theory," 
i.e.  the  expressions  are  generalizations  of  those  conventionally 
used  for  straight  beams  when  shear  effects  and  stresses  at 
right  angles  to  the  axis  of  the  beam  are  ignored.  The  approach 
is  similar  to  the  treatment  of  straight  beams  of  thin-walled 
section  in  Ref.  1,  Chapter  5,  where  the  sections  are  considered 
to  be  built  up  from  a number  of  flat  plates  for  each  of  which 
the  strain  energy  from  elementary  beam  theory  is  known.  Making 
the  assumption  that  the  shape  of  the  cross  section  does  not 
change,  the  energy  in  each  plate  can  be  expressed  by  the  global 
coordinates  of  the  bar.  Allowing  for  the  continuity  of  strains 
wherever  two  of  the  plates  are  jointed,  strain  energy  expressions 
and  the  location  of  the  shear  center  are  obtained. 

In  the  present  treatment  the  elements  of  which  the  bars 
consist  are  either  not  flat,  or  not  straight.  The  sections  shown 
in  Figs,  la,  b,  and  c will  be  treated.  The  ring  in  Fig.  la  is 
of  simple  rectangular  cross  section;  it  is  really  a flat  plate 
of  annular  shape.  In  Fig.  lb  the  section  is  a T.  The  web  is 
again  an  annular  flat  plate,  while  the  flange  is  a short  segment 
of  a cylindrical  shell.  The  third  case,  Fig.  lc,  consists  of 
three  parts  of  similar  nature.  Expressions  for  the  strain 
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energies  of  the  elements  are  derived  in  the  Appendices,  using 
appropriately  simplified  relations  available  from  plate  and 
shell  theories.  For  deformation  of  the  entire  bar  in  the 
plane  of  curvature  a generalization  of  conventional  beam  theory 
is  used,  which  assumes  that  plane  sections  remain  plane,  and 
also  at  right  angles  to  the  deformed  center  line. 

The  expressions  in  the  Appendices  could  also  be  used  to 
treat  nonsymmetric  cases  like  U-  or  L-  stiffeners.  Figs.  Id,  e. 
These  cases  are  not  included  because  their  use  seems  rather 
unlikely. 

The  strain  energy  expressions  V obtained  may  be  used 
in  various,  fairly  obvious  ways.  They  may  be  utilized  to 
find  the  boundary  conditions  at  the  stiffeners  for  the  a priori 
known  partial  differential  equations  for  vibrations  of  the  shell 
by  using  Hamilton's  principle, 

/(v  - T) dt  = extremum  (1) 

and  applying  calculus  of  variation  with  respect  to  the  shell 
coordinates  x and  <J>,  Fig.  2.  V and  T are  the  strain  energy 
and  the  kinetic  energy,  respectively. 

As  an  alternative,  one  can  use  a Raleigh-Ritz  approach 
and  introduce  appropriate  approximation  for  the  shell  displace- 
ments u,  v,  w.  Fig.  3,  into  Eq.  (1). 

As  a further  alternative,  one  may  introduce  into  Eq.  (1) 


the  expressions 
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u = U (x)  cos  (ni  + a) 
n 

v = Vn(x)  sin(n<j>  + a)  (2) 

w = W (x)  cos(n4>  + o) 
n 

where  a is  a phase  angle.  This  substitution  reduces  the  partial 
differential  equations  of  the  shell  to  three  ordinary,  simul- 
taneous ones  in  U,  V and  W. 

For  either  of  the  alternatives,  a suitable  expression 
for  the  strain  energy  of  the  cylindrical  shell  may  be  found  in 
Ref.  2. 
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II.  STIFFENING  RINGS  OF  RECTANGULAR  CROSS  SECTION,  FIG.  la 

Fig.  4a  shows  a portion  of  the  shell  of  thickness  t and 
radius  a and  an  interior  stiffener  of  depth  d and  thickness  h. 
The  depth  d is  a nominal  one,  measured  from  the  center  surface 
of  the  shell  to  the  innermost  edge  of  the  stiffener.  Fig.  4a. 
This  figure  also  shows  the  centroid  0 of  the  stiffener  cross 
section,  and  the  radius  RQ  of  the  centroidal  circle.  Fig.  4b 
shows  the  original  and  the  displaced  center  lines  of  shell 
and  stiffener,  and  the  displacements  uQ,  wQ  of  the  centroid  0 
as  well  as  the  rotation  8.  The  out-of-plane  displacement  is 
vQ,  but  cannot  be  indicated  in  Fig.  4b. 

Using  Eqs.  (A-13)  and  (A-22)  for  the  portions  of  the 
strain  energy  of  the  stiffener  in  and  out  of  the  plane  of 
curvature,  respectively. 


1 F 7 32wn  2 FA  9V0  2 

v-?ni  v + w°’  +!(u°  + ^’  + 


(3) 
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where 


7 _ hd3  . _ . . T _ h3d  T _ h3d  D _ a d ... 

Z 12  ' A hd'  Io  12  ' J0  3 ' R0  2 (4) 


The  value  of  Z is  approximate,  but  suitable  if  d <<  a.  The 
displacements  uQ,  w^  and  8 can  be  expressed  from  geometry  by 
the  shell  displacements  at  point  A,  Figs.  4a,  4b, 
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U0  = UA  - 2 8 = UA  + 2 -3F  (5) 


8vo 

The  quantity  vQ  appears  in  Eq.  (3)  only  in  the  form  . 

To  express  this  derivative,  the  equality  of  the  strains  in  shell 
and  stiffener  at  point  A,  Fig.  4a,  is  used 
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The  left  hand  side  of  this  equation  is  the  membrane  strain  in 
the  shell  at  point  A,  while  the  right  hand  side  is  obtained 
from  Eq.  (A-15)  for  g = d/2.  Noting  the  relation  between  a and  RQ 
gives 


9vo  _ d ®2"a  ro  9va 

~ 2a~^  a a* 


(7) 


Substitution  into  Eq.  (3)  gives  the  strain  energy  in  the  ring 
stiffener  in  terms  of  uA,  v^  and  wft. 
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If  desired,  one  could  introduce  the  approximation  RQ  a,  but 
the  resulting  simplification  is  hardly  worthwhile. 
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III.  STIFFENING  RINGS  OF  T-SECTION,  FIG.  lb 

Figure  5a  shows  the  shell  and  an  interior  stiffener 
and  all  dimensions.  The  centroid  of  the  stiffener  is  at  a 
distance  e from  the  middle  surface  of  the  shell.  Figures  5b 
and  5c  show  the  web  and  flange  in  their  original  and  in 
their  displaced  positions,  respectively.  Also  shown  are  the 
centroids  C>w  and  01  and  their  displacements. 

The  strain  energy  of  the  entire  stiffener  due  to  dis- 
placements wQ  and  vQ  is 


V(v0'W0)  = I 


]d * 


(9) 


where  ZQ  = I and  AQ  are  the  moment  of  inertia  and  the  area  of 
the  section,  respectively. 

The  strain  energies  due  to  the  displacement  in  the 
x-direction  and  rotation  for  the  web,  V , and  for  the  flange, 

, are  according  to  Eqs.  (A-22)  and  (A-32) , respectively, 
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where  the  section  properties  are  defined  by 


I = 
w 


dh' 


12 


J = 
w 


dh' 


X1  " 


fclb' 

12 


j 

J1  " 3 


(12) 


The  displacements  wQ,  and  uw  can  be  expressed  by  the 
equivalent  quantities  at  point  A.  The  relations  are 


d 9WA 

”o  - "s'  uw  * “a  - 2 8-  U1  = UA  - ds-  e ' - -5F 


(13) 


8vo  . 

To  express  . . , the  strain  in  the  shell  at  point  A,  and  the  strain 

o (p 


in  the  web  at  the  same  point  are  equated 

1 1 ^ V A 1 ^ v0  e 

a WA  + a T$  = (W0  + Hh0  ~ R0(RQ+e) 


32w, 


3 4> 


i + V 


(14) 


The  right-hand  side  of  this  equation  is  Eq.  (A-15)  for  n = e. 
After  simplification 
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The  total  strain  energy  is  the  sum  of  Eqs.  (9) , (10)  and 
(11).  After  substitution  of  Eqs.  (13)  and  (15) 
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This  expression  can  be  somewhat  simplified,  if  desired, 
by  noting  that  d <<<  a,  and  that  for  thin-walled  sections 
Iw  <<<  1 2.*  Using  also 


R. 


0 — Ri  — Rw^.aZc3^.a-2d:~a 


one  obtains 


2 

1 3 wa  o 3 w.  3v.  2 

7 - -7  ;IEI0<— f+  V + eV“'a  + « — I +*—  > + 

1 34>  94)^  94> 


2a' 


♦EI^ 


92qa  93w.  9w  _ 

d ^-a^-)2+  G(J..+J,)(^+a 


9u, 


.2 
9 w 


(17) 


A,  2 


94>  9x94>‘ 


w ' 94,  - 


1 CUP 


-9- 


IV.  I-BEAM  STIFFENING  RINGS,  FIG.  Ic 
Figure  6a  shows  a cross  section  of  the  shell  and  of 
the  stiffening  ring  inside  the  shell.  The  centroid  of  the  web 
and  of  the  entire  stiffener  is  at  0.  The  properties  of  the 
web  and  its  displacements  will  carry  the  subscript  0.  The 
two  flanges  are  of  equal  dimensions,  b by  t^,  and  their 
centroids  are  0^,C>2,  respectively.  The  displacements  of  the 
two  centroids  carry  the  subscripts  1,  2,  respectively. 

Just  as  in  Sections  I and  II,  the  portion  of  the  strain 
energy  V(wQ,v0)  is  given  by  Eq.  (A-13) , 


1 r rEZ  , 3 W0 

2 f[~3  ( 2 

Ro  9<)) 


EA  , 

^ (w° 


where  A and  Z are  the  area,  and  the  moment  of  inertia  of  the 
entire  section. 

The  portion  VQ  of  the  strain  energy  of  the  web  due  to 
the  displacements  uQ  and  6 is  according  to  Eq.  (A-22) 


o 2 

7 + V>  + 


G J 3 u 

0,  0 

*3  1 94> 
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where 


(20) 


The  portions  of  the  strain  energy  due  to  the  displacements 
of  the  two  flanges  u^  and  6 are,  respectively,  from  Eq.  (A-32) 
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V1(u1,6) 
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where  u^,u2  are  the  displacements  of  the  flanges  1 and  2 
respectively,  and 


I = 


tlb~ 

12 


J = 


fclb 


R1=a-e,  R2  = a - e - d (23) 


Referring  to  Fig.  6b,  the  displacements,  except  , 
can  be  expressed  by  the  displacements  of  the  shell  at  point  A, 


9wa 
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u2  = uA  -(d  + e)6=  uA+  (d  + e) 
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The  equality  of  the  hoop  strains  at  A in  the  shell  and 
in  the  stiffener,  with  n = e + d/2  gives 


d + ze  9 w 
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and  after  simplification 
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Using  Eqs.  (24)  and  (25),  the  four  Eqs.  (18)  , (19)  , (21) 


and  (22)  become 
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(26) 


The  entire  strain  energy  V is  the  sum  of  the  four  terms, 

Eqs.  (26) . However,  some  approximations  seem  permissible.  The 
first  term  in  the  relation  for  contains  the  factor  IQ  = dh^/12, 

which  is  very  small  in  comparison  to  the  similar  terms  in  and  , 

where  I = t^b3/12.  For  thin-walled  sections  the  term  multiplied 
by  IQ  is  negligible.  Further  e = (t  + t^)/2  is  small  compared  to 
d or  R.  Using  a ~ 0 is  thus  a reasonable  approximation.  Assuming 
further  that  d <<<  a,  one  may  use  approximately  R0  ^ Rj  ; R2  i a* 

The  end  result  is 
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V.  SUMMARY 

Strain  energy  expressions  have  been  obtained  for  ring 
stiffeners  of  the  types  shown  in  Figs,  la,  b and  c.  The  results 
are,  respectively,  given  by  Eqs.  (8),  (17)  and  (27).  The 
stiffeners  shown  in  Fig.  1 are  all  on  the  inside  of  the  shell. 
However  the  results  obtained  are  also  applicable  to  stiffeners 
of  the  same  type  on  the  outside  of  the  shell.  In  such  cases  the 
quantities  d and  e occurring  in  Eqs.  (8),  (17)  and  (27)  must  be 
replaced  by  -d,-e,  respectively. 

Possible  applications  for  the  strain  energy  expressions 
are  indicated  in  the  Introduction.  The  suggested  uses  require 
an  expression  for  the  strain  energy  of  shell  panels  adjoining 
the  stiffeners.  Such  an  expression  is  available  in  Ref.  2. 
Defining  the  shell  displacements  u,v  and  w as  shown  in  Fig.  3, 
the  strain  energy  of  a panel  is 


^ ff[ a2u2  + (v^  + w)  2 + 2avux(v^  + w)  + ^^(u^  + avx)2]dxd<J> 


— — 2 + (w  + w)  2 + —5— -(  (aw  + u )2 

n .*s  3 XX  4®  2 xd>  <t> 


24(1-0  a‘ 


2 2 

(v  “ w„a)  + 2va  w (w. 
x xq>  xx  d 


- v. ) - 2a  u w ]dxd(j) 
4>  x xx 


where  t is  the  shell  thickness,  and  the  subscripts  x or  4>  indicate 
Partial  derivatives  with  respect  to  x or  *)>.  The  double  integrals 
extend  over  the  area  of  the  shell  panel. 
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APPENDIX  1. 

Strain  energy  of  a circular  curved  bar  due  to  deformation  v,  w 
in  the  plane  of  curvature.  Fig.  A-l. 

Consider  an  element  of  a curved  bar.  Fig.  A-2,  the  cross 
section  of  which  is  symmetric  to  the  plane  of  curvature.  The 
element  is  stressed  by  a moment  M and  a direct  force  N,  both 
in  the  plane  of  curvature.  The  small  deflections  of  a point  P 
in  the  location  <J>  can  be  described  by  the  radial  component  w and 
the  tangential  component  v.  Fig.  A-l.  Figs.  A-l  and  2 show  I-beams 
with  unequal  flanges,  which  are  generalizations  of  the  three 
specific  cases.  Figs,  la,  b and  c,  needed  in  the  body  of  the 
report. 

The  approach  used  is  a generalization  of  the  simple 
Navier-Euler  theory  in  straight  beams,  which  assumes  that  only 
stresses  in  the  axial  direction  contribute  materially  to  the 
strain  energy.  Stresses  in  the  radial  direction  and  shear 
stresses,  both  of  which  must  exist,  are  thus  assumed  to  con- 
tribute only  negligibly  to  the  strain  energy.  The  assumption 
that  plane  cross  sections  remain  plane  and  at  right  angles  to 
the  deformed  center  line  of  the  bar,  leads  to  the  distribution 
of  bending  stresses  as  obtained  in  the  classical  treatment  of 
Winkler- Resal.  The  following  is  not  concerned  with  a 
re-derivation  of  the  equations  for  the  stresses,  but  with  the 
formulation  of  strain  energy  expressions  in  terms  of  derivatives 
of  v and  w.  Such  expressions  seem  not  to  be  available  in  the 


literature. 
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Consider  an  element  of  the  bar  of  length  ds  in  its 
original  and  in  its  distorted  shape.  Fig.  A-2.  The  centroidal 
axis  of  original  length  ds  will  be  lengthened  by  Ads,  and  the 
angle  between  the  two  faces  will  change,  as  shown  in  the  same 
figure,  by  + Ad<J>. 

Considering  an  element  dA  at  a radial  distance  n from 
the  centroid,  one  can  compute  strain  and  stress  in  the  element 
as  function  of  Ads  and  Ad<J> 


o 


P,1  Ads  n Ad<fr  . 
d<j>  R + n d<J>  ' 


(A-l) 


As  usual,  the  quantity  Z is  introduced 


Z 


R / 


r + n 


dA 


r2  r ndA 
R J r + n 


(A-2) 


If  the  depth  of  the  bar  is  small  compared  to  the  radius  R,  the 
value  of  Z is  practically  identical  to  the  moment  of  inertia  I. 
In  the  body  of  the  report,  it  will  be  assumed  that  Z I. 
Comparing  the  resultants  of  the  stresses  given  by  Eq.  (A-l) , 
and  M and  N,  one  can  determine  Ads  and  Ad <J> , 


Ads  _ R 
“3$  EA 


(A- 3) 


The  strain  energy  dV  in  the  element  of  length  ds  = Rd4> 
being  equal  to  the  work  done  by  the  forces  M and  N during  the 
distortion  of  the  element,  one  finds 
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V = /dV  Rd»  = i /|f  <^)2  + f (^)2|d*  (A-4) 

The  quantities  and  are  to  be  expressed  in  terms 

of  v and  w,  which  are  the  components  of  the  displacement  of  the 
centroid  0 of  the  cross  section.  Figuare  A-3  shows  the  original 
and  the  distorted  element  superimposed  on  each  other.  Using 
polar  coordinates,  p(<p)  = R + w(<J>),  the  curvature  of  the  original 
center  line  is  1/R,  while  the  curvature  1/R^  of  the  distorted 
center  line  is,  with  p = R + w. 


p2+2,|£,2  -Pa!| 

1 sLfc_ 


9 2 3/2 

ip2*‘fs*  I 


1 

R+w 


1 A2 

1 d w 

2 2 
( R+w ) * d * 


(A-5) 


The  approximate  result  is  obtained  by  using  the  fact  that 
dw  2 2 

<<  (R+w)  . Forming  the  expression  1/R^^  - 1/R,  and  allowing 
for  w <<  R,  and  ds  = Rd4>,  one  finds 


1 

R 


w 

R(R+w) 


1 

(R+w)2 


+ w) 


(A- 6) 


In  addition  to  above  relation  there  is  a geometric  one 
between  R,  R^ , Ads  and  Ad<J>,  which  can  be  read  from  Fig.  A-3. 

The  total  length  of  the  distorted  axis,  ds  + Ads,  must  equal  the 
new  radius  Rj  multiplied  by  the  angle  a enclosed  by  the  two 
faces  of  the  deformed  element.  Thus 
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Rla  R1(X  + Ad*  + = As  + Ads  (A-7) 

Rearranging  and  dividing  by  R^R  d(()  gives 

1 1 _ 1 Ad(J)  Ads  R1~R  1 Ad£  (A-8) 

R1  ~ R R d4>  + "3T  r2r_  - R d4> 

The  approximation  in  Eq.  (A-8)  is  permissible  because  (R^-R)/R 
is  inherently  a small  quantity  in  comparison  to  unity.  Equations 
(A-6)  and  (A-8)  furnish 


1 ,d2w 

R 57 


+ w) 


( A-9 ) 


In  conjunction  with  Eq.  (A-3)  this  relation  leads  to  the  well-known 
differential  equation  for  the  radial  displacement  w,  see  Ref.  (3) . 

An  additional  geometric  relation  can  be  obtained  from 
Fig.  A-4.  The  end  points  A and  B of  the  element  ds  displace 
to  A'  and  B',  respectively.  The  distances  CA'  and  CB ' follow 
from  Fig.  A-4,  where  quantities  which  are  small  of  higher  order 
are  neglected. 


CA1  = dw  - vdiji  - dvd<))  ^ dw  - ^ vds 

( A-10 ) 

CB  * = Rd$  + wd<)>  + dv  + dwd<{>  ~ ds  + ^ wds  + dv 

R 


Further 
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dw 

ds 
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Expanding  the  square  root  by  the  binomial  law  and  neglecting 
higher  order  terms  gives 


Ads  . dv  . . 

-a?  = ” + a*  (A-12) 

Substitution  into  Eq.  (A-4)  gives  finally  the  strain  energy 
expression 

2 2 2 
1 ,, EZ  ,d  w , . . EA  . , dv.  , 

V = ~ j + w)  + — (w  + — ) ] d<t»  (A- 13) 

R d«j>  R d<(t 

It  is  noted  that  the  second  Eq.  (A-3)  and  Eq.  (A-12) 
give  the  differential  equation 


w dv 
R ds 


EA 


IN  + §] 


(A-14) 


Conventional  texts  contain  only  an  approximation  of  the  equation 
where  the  term  M/R  does  not  appear. 

It  will  also  be  necessary  to  have  an  expression  for  the 
strain  in  a location  n,  Fig.  A-2.  Using  Eqs.  (A-l) , (A-9)  and 
(A-12)  one  finds 
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1 Ads  n Ads  _ 1_ . dv. 

R d4>  R+n  d4>  RtW  d<J> 
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r ( R+n ) 


,2 

,d  w 

' 2 
d<p 


+ w) 


(A-15) 
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APPENDIX  2. 

Strain  energy  of  a thin  annular  plate,  displaced  at  right  angles 

to  its  middle  plane.  Fig.  5. 

The  web  of  a ring  stiffener  of  a shell.  Fig.  la,  b or  c, 
may  be  considered  as  an  annular  plate.  An  expression  for  the 
strain  energy  of  such  a plate  in  polar  coordinates  can  be  found 
in  Ref.  (4  p.  346,  Eq.  (o) ) . This  expression  is  a double 
integral  over  three  major  terms.  The  sum  of  the  first  two  terms 
can  be  recognized  as  due  to  the  direct  stresses,  while  the  third 
is  due  to  the  shear  stresses.  Let  V = + V2 , and  using  the 

symbol  u for  the  normal  displacements  of  the  plate. 


V,  = 


? ~ 
2 3r2  r 9r 


i *2-  2 

1 3 u. 

1 — 2 ~ 

r 34> 


-2  ( 1-v) 


.2- 
3 u 


3r 


(I  3u  + l_ 

r 8r  r2 


a2-  2 

— y)  ] rdrd4> 
34> 


(A-16a) 


v)// 


,1  32u  1 32G,2 

^ ' 7 17 


rdrd4> 


(A-16b) 


Consistent  with  the  assumption  that  the  only  stresses 
contributing  materially  to  the  strain  energy  of  bending  are 
direct  stresses  a.,  the  term  D in  (A-16a)  is  to  be  evaluated 
for v = 0,  while  the  value  v / O is  retained  in  Eq.  (16-b)  for 
shear  effects.  Further,  as  in  elementary  beam  theory,  it  is 
assumed  that  the  radial  axis  of  the  cross  section  remains  a 
straight  line 
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u (r)  = u (n)  = u + n6 


(A-17) 


where  u is  the  displacement  of  the  centroid  and  6 the  rotation 
of  the  cross  section.  Using  r = R + n # dr  = dn,  and  Eq.  (A-17), 
one  has  the  relations 


3u  _ 3u,  _ R 
3r  - 3n  " b 


.2 
9 u 

9((>2 


+ n 


3u  _ 3u  3^u  _ j|J3 

d<P  3<f>  ^ 3<j>'  9r34»  3<t> 


(A-18) 


After  substitution  of  these  relations  into  Eqs.  (A-16) 
one  can  integrate  with  respect  to  n = r - R.  Using  the  fact 
that  d/ R is  small  compared  to  unity,  the  various  integrals 
can  be  approximated 
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r dr 
f ~ 

r 


= £n 
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' 1 ~T 

r 


2 ( R-d/2 ) 


2 2 
R -d^/4 


d 

^7 


(A-19) 


One  obtains  thus 


1 Edh" 

1 - 2 i2 
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and 


V, 


1 Gdh3  ,1  , 3u 

2 3 1 ^ l30 


p ag.2 

R 


( A- 2 1 ) 


As  a check  on  the  adequacy  of  the  simplifications  made 
in  Eqs.  (19)  consider  a rigid  body  rotation  of  the  annular  plate 
with  respect  to  the  axis  0 = + t/2,  0 = - tt/2,  where  u = A cos  <J> , 

8 = ^ cos  0 and  A is  the  displacement  at  0 = 0.  Substitution 
into  Eqs.  (A-20,  21)  indicates  that  for  this  displacement 
= V 2 = as  required  for  a rigid  body  displacement. 

The  total  strain  energy  for  out-of-plane  displacements 

is  thus 


V = 


Gdh' 


(S  - R S ld0  (A- 22) 


38 


30 


¥0 


The  two  coefficients  appearing  in  the  integrand  are  the  values 
El  and  GJ  appearing  in  the  equivalent  expressions  for  a straight 
beam.  For  R -*■  <*> , Eq.  (A-22)  thus  furnishes  the  conventional 
value  for  V for  a straight  beam. 
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APPENDIX  3. 

Strain  energy  stored  in  the  flanges  of  rings  of  T-  and  I-Sections 

Let  the  displacements  of  an  arbitrary  point  P on  the  middle 
surface  of  a flange,  Fig.  7,  be  designated  by  the  symbols 
UP'VP,WP  the  displacements  of  the  centroid  0 of  the 

undeformable , rectangular  cross  section  are  u,v,w  and  the  rota- 
tion is  6. 

Only  sections  which  are  symmetric  to  the  plane  of  curvature 
are  treated  so  that  the  strain  energy  separates  into  the  sum  of 
two  terms 


V = V1 (w,v)  + V2 (u,B) 


(A-23) 


one  term  depending  only  on  w and  v,  the  other  on  u and  3.  The 
term  has  already  been  obtained  in  Appendix  1,  Eq.  (A-13) 
and  the  associated  strain  in  Eq.  (A-15) . 

The  second  term,  V2 , for  T-  or  I-sections  will  be  derived 
in  the  body  of  the  report,  using  an  expression  for  the  web  alone 
derived  in  Appendix  2,  in  conjunction  with  a relation  to  be 
derived  here,  treating  flanges  of  rectangular  cross  section, 
b x t,  as  short  pieces  of  cylindrical  shells.  See  Fig.  A-7. 
Assuming  that  the  cross  section  does  not  change  its  shape, 
the  displacements  of  a point  P in  the  location  n , z = 0,  Fig.  A-7, 
are 


w 


P 


nB,  Up  = u vp 


n 

R 3<f> 


( A-24 ) 
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The  strain  energy  V2(u,0)  can  be  divided  in  a portion 

due  to  hoop  strains  e,.,  and  one  due  to  shear  strains 

y ,V2  = V2t  +V2y*  Allowing  again  for  the  fact  that  direct  stresses 

other  than  og  are  negligible,  the  value  of  Poisson's  ratio  in 

V_  is  assumed  to  vanish 
2e 

V2f  = | fff  dndz(R  + z ) d<f>  (A-25) 

while  the  usual  value  of  v is  retained  in 

^2y  = J Iff  y2dndz(R  + z)d<j>  (A-26) 

When  evaluating  Eq.  (A-25)  it  is  assumed  that  the  bending 
strains  do  not  vary  significantly  through  the  thickness,  t <<  d, 
t <<<  R,  so  that 


V2e  = T'  //eS  dn(R  + z)d * - //eS  dnd4)  (A-27) 

where  the  value  cc  at  z = 0 is  to  be  used.  Reference  [5,  Eq.(5-b) 
on  p.  209]  gives  for  z = O,  after  substitution  of  Eqs.  (A-24) 


. 3v  w 
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S R 3<J>  R 
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R 34>  R 


Equation  (A-27)  gives  thus  finally 
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(A-28) 
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The  above  integral  vanishes,  as  required,  for  the  rigid  body 

motion  u = A cos  <J> , 8 = 4 cos  <J>. 

To  evaluate  the  integral  in  Eq.  (A-26)  use  is  made  of 

Ref.  (5,  Eq.  (5-c)].  The  value  of  the  shear  strain  at  an  arbitrary 

point  A,  Fig.  A-7,  is  expressed  by  the  values  at  points  P on  the 
center  plane 


■y 
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R+z 


9u 


R+z 


R 9<J> 


9x9  <p 


(-  + — ) 
R+z 


( A-29 ) 


Substitution  of  Eqs.  (A-24)  gives 


2 R+z  9u 
R2  (R+z)  d(f> 


z 2R+Z  iB  _ t z f 1 9u  _ 1 9B 
R(R+z)  9<,)  ~ R2  9<|>  R 9 <y 


( A-30) 


The  approximation  used  utilizes  the  fact  that  max  z = t/2  <<<  a. 
The  value  of  the  integral  in  Eg.  (A-26)  becomes  thus 


1 Gbt3 

2 3 


d4> 


(A- 31) 


As  necessary,  this  value  vanishes  for  the  rigid  body  motion 
tested  on  Eq.  (A-28) . 

The  total  value  of  V2(u,B)  is 


V2  (u,B] 


1 Etd3  r 1 92u  . do v 2 , . , 1 Gbt3  .1  . 9u  0 98  2 . 

2 1 — 2 + R8)  d<t>  + 7 f ■“ T (aX  ' R TT) 

12  RJ  23  R3  9<t) 


(A-32) 


In  the  limit,  R -*■  <»,  Eq.  (A-32)  furnishes  the  usual 
expression  for  the  strain  energy  of  a straight  bar. 
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} 13  ABSTRACT 


Strain  energy  expressions  are  obtained  for  rings  of 


’ rectangular,  T-  and  1-  section.  The  expressions  are  intended  for 
1 

> use  in  the  dynamic  analysis  of  ring  stiffened  cylindrical  shells.. 

I 

* The  approach  is  essentially  a generalization  of  the  conventional, 

l 

approximate  analysis  of  straight  beams,  i.e.  the  influence  of 

p 

shear  stresses,  and  of  direct  stresses  at  right  angles  to  the  axis 
. of  the  beams  is  neglected. 
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